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ABSTRACT 

We present measurements of the intergalactic medium (IGM) temperature within 
^ 5 proper Mpc of seven luminous quasars at z ~ 6. The constraints are obtained 
from the Doppler widths of Lya absorption lines in the quasar near-zones and build 
upon our previous measurement for the z — 6.02 quasar SDSS J0818-I-1722. The 
expanded data set, combined with an improved treatment of systematic uncertainties, 
yields an average temperature at the mean density of log(To/K) = 4.21±g Q3 (±g Qy) 
at 68 (95) per cent confidence for a flat prior distribution over 3.2 < log(ro/K) < 4.8. 
In comparison, temperatures measured from the general IGM at z ~ 5 are ~ 0.3 
dex cooler, implying an additional source of heating around these quasars which is 
not yet present in the general IGM at slightly lower redshift. This heating is most 
likely due to the recent reionisation of Hen in vicinity of these quasars, which have 
hard, non-thermal ionising spectra. The elevated temperatures may therefore represent 
evidence for the earliest stages of He ii reionisation in the most biased regions of the 
high-redshift Universe. The temperature as a function of distance from the quasars is 
consistent with being constant, log(To/K) ~ 4.2, with no evidence for a line-of-sight 
thermal proximity effect. However, the limited extent of the quasar near-zones prevents 
the detection of Hem regions larger than ~ 5 proper Mpc. Under the assumption the 
quasars have reionised the He ii in their vicinity, we infer the data are consistent with 
an average optically bright phase of duration in excess of lO^ '^ yr. These measurements 
represent the highest-redshift IGM temperature constraints to date, and thus provide 
a valuable data set for confronting models of H i reionisation. 
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1 INTRODUCTION 

The epoch of reionisation marks the time when luminous 
sources became a dominant influence on the physical state 
of the intergalactic medium (IGM). Studying the impact 
of reionisation on the IGM therefore provides valuable in- 
sight into the formation and evolution of the flrst stars and 
galaxies (see Morales & Wyithe 2010 for a recent review). 
One quantity which is closely related to the timing and 
extent of the reionisation epoch is the temperature of the 
low-density IGM. The IGM temperature is set primarily by 
photo-heating and adiabatic cooling at z < 10, and therefore 
depends on the spectral shape of the radiation from early 
ionising sources and the time elapsed since intergalactic hy- 



drogen and helium were reionised [e.g. Miralda-Escude & 
Rees 1994; Theuns et al. 2002a; Hui & Haiman 2003). 

Existing measurements of the IGM temperature are ob- 
tained from the forest of intergalactic Hi Lya absorption 
lines observed in the spectra of bright quasars. The widths of 
the Lya absorption lines are sensitive to the temperature of 
the IGM through a combination of thermal (Doppler) broad- 
ening and pressure (Jeans) smoothing of the underlying gas 
distribution [e.g. Haehnelt & Steinmetz 1998; Peeples et al. 
2010), in addition to broadening from peculiar motions and 
the Hubble flow [e.g. Hernquist et al. 1996; Theuns et al. 
2000). Consequently, using a statistic sensitive to the ther- 
mal broadening kernel combined with an accurate model for 
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measurement calibration (typically high resolution hydro- 
dynamical simulations of the IGM), various authors have 
placed constraints on the thermal evolution of the IGM at 
2 < z < 4.8 using the Lya forest (Schaye et al. 2000; Ricotti 
et al. 2000; Thcuns & Zaroubi 2000; McDonald et al. 2001; 
Zaldarriaga 2002; Lidz et al. 2010; Becker et al. 2011). 

Becker et al. (2011) recently presented the most precise 
IGM temperature measurements to date obtained from the 
Lya forest using a set of 61 high-resolution quasar spec- 
tra over a wide redshift range, 2 < z < 4.8. These au- 
thors found that the IGM temperature at mean density, 
To, is ~8000K at z ~ 4.4 and gradually increases toward 
lower redshift. Becker et al. (2011) attributed this heat- 
ing to the onset of an extended epoch of Hcii rcionisation 
driven by (juasars, which produce the hard ionising photons 
necessary for doubly ionising intergalactic helium (Madau 
et al. 1999; Purlanetto & Oh 2008; McQuinn et al. 2009). 
When combined with measurements of the intergalactic 
Hen Lya opacity at 2<3 (e.g. Fechncr et al. 2006; Shull 
et al. 2010; Worseck et al. 2011; Syphers et al. 2011), these 
data suggest Hen reionisation was in its final stages by or 
after 2 ~ 3. 

However, in order to probe deeper into the epoch of 
H I reionisation at 2: > 6, constraints on the IGM tem- 
perature at higher redshift are required. Unfortunately, 
due to the increasing Lya opacity of the IGM (Songaila 
2004; Fan et al. 2006; Becker et al. 2007), measurements 
of the IGM temperature from the general Lya forest are 
extremely challenging at 2; > 5. Bolton et al. (2010) (here- 
after BIO) recently side-stepped this difficulty by measur- 
ing the IGM temperature in the vicinity of the z = 6.02 
quasar SDSS J0818-I-1722. Using detailed numerical simu- 
lations, BIO demonstrated that the cumulative probability 
distribution function (CPDF) of Doppler parameters mea^ 
sured from Lya absorption lines in the highly ionised near- 
zone is sensitive to in-situ photo-heating by the quasar, as 
well as heating from earlier ionising sources. BIO therefore 
used the Doppler parameter CPDF to infer the IGM tem- 
perature at mean density within ~ 5 proper Mpc of the 
quasar, log(To/K) = 4.371q;°5, P®"" ^^^^^ confidence. 

However, the analysis of only a single linc-of-sight means 
the constraint presented by BIO has a rather large statisti- 
cal uncertainty, and is not a reliable representation of the 
IGM thermal state at z ~ 6 as a whole. Temperature mea^ 
surements from independent lines-of-sight would therefore 
significantly aid in reducing the statistical uncertainty on 
an averaged measurement, as well as enabling an investiga- 
tion of line-of-sight temperature variations. 

The goal of this study is to extend and improve upon the 
preliminary work of BIO in several important ways. Firstly, 
we analyse the line widths in the near-zones of six further 
quasars at a ~ 6 using high-resolution data obtained with 
Keck/HIRES and Magellan/MIKE (Becker et al. 2006, 2007; 
Calverlcy et al. 2011). These arc combined with improved 
measurements of the quasar systemic rcdshifts and absolute 
magnitudes (Carilli et al. 2010). Secondly, we significantly 
expand the suite of numerical simulations used to calibrate 
our temperature measurements, enabling us to explore a 
wider range of thermal histories. Finally, we undertake a 
more detailed analysis of the key systematic uncertainties 
identified in BIO: metal line contamination and continuum 
placement. In addition, we also now investigate the impact 



of the uncertain thermal history of the IGM at z > 6 on the 
temperature measurements. 

This paper is structured as follows. In Section 2 we 
present the observational data and numerical simulations 
used in this work. In Section 3, we briefly review the method 
for obtaining the near-zone temperature constraints, which 
was originally discussed in detail by BIO. We investigate the 
important systematic uncertainties on our measurements 
in Section 4, and in Section 5 we present our tempera- 
ture measurements before finally concluding in Section 6. 
An appendix which presents some tests of our temperature 
measurement procedure is given at the end of the paper. 
Throughout we assume fim = 0.26, Oa = 0.74, Sib = 0.0444, 
h = 0.72, CTg = 0.80, ris = 0.96 (Komatsu et al. 2011) and 
a helium fraction by mass of F = 0.24 (Olive & Skillman 
2004). All distances are quoted in proper units unless oth- 
erwise stated. 



2 DATA AND NUMERICAL MODELLING 
2.1 Observational data 

We analyse seven high-resolution quasar spectra in this 
work. The spectra for six of these objects were obtained us- 
ing the High Resolution Echelle Spectrograph (HIRES, Vogt 
et al. 1994) on the 10-m Keck I telescope using a 0.86"slit, 
which produces a resolution oi R = 40 000 (FWHM = 
6.7kms~^). The seventh spectrum was obtained using the 
Magellan Inamori Kyocera Echelle (MIKE, Bernstein et al. 
2003) spectrograph on the 6.5-m Magellan II telescope. 
This instrument has a slightly lower resolution compared to 
HIRES, with FWHM=13.6 kms"^ The data were processed 
using a custom set of IDL routines (Becker et al. 2006, 2007) 
that include optimal sky subtraction (Kelson 2003). The fl- 
nal binned pixel sizes are 2.1 km s^^ for the HIRES spectra, 
and 5.0 km s"^ for the MIKE data. 

The quasars are fully summarised in Table 1, and the 
spectra are displayed in Figures 1 and 2 . The systemic red- 
shifts and absolute AB magnitudes are taken from Carilli 
et al. (2010). The continuum fitting in the near- zones was 
performed by first dividing the spectra by a power-law with 
fv oc 1/"° ®, normalised near 1280 (1 -|- z) A. The results of 
this procedure are displayed in Figure 1. The Lya emission 
line was then fitted with a slowly varying spline, shown by 
the dotted red curves in Figure 1. The fully normalised near- 
zone spectra are displayed in Figure 2. Further details on the 
continuum fitting procedure may be found in BIO, and the 
impact of the continuum fitting uncertainties on our results 
is examined in detail in Section 4.2. 

We obtain the Doppler parameter CPDF for the 
Lya absorption lines following a procedure identical to the 
approach described by BIO. We use an automated version 
of the Voigt profile fitting package VPFIT^, to fit lines 
over the velocity ranges indicated in Table 1. All lines with 
log(A^Hi/cm~^) > 17, 6 > lOOkms"^ and with relative er- 
rors in excess of 50 per cent are discarded for the tempera- 
ture analysis. These discarded lines make up 30 per cent of 
the total fitted to the data. The remaining line fits are indi- 
cated by the red dotted curves in Figure 2, with the centre of 

^ http://www.ast.cam.ac.uk/~rfc/vpfit.html 
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Figure 1. The seven quasar spectra analysed in this work against velocity with respect to the rest frame Lya wavelength. The spectra have 
already been normalised by a power-law with fi, oc u~'^'^ (see text for details). A slowly varying spline (dotted red curves) has also been 
fitted to the Lya emission lines in each spectrum to approximate the quasar continuum. A detailed analysis of the systematic uncertainty 
associated with the continuum fitting procedure is presented in Section 4.2. 
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Figure 2. Transmitted fraction against tlie Hubble velocity blueward of the Lya emission line for the seven quasar near-zones analysed 
in this work. The red dotted curves display the Voigt profile fits to the absorption lines made using VPFIT, with the centre of the lines 
marked by the downward pointing arrows. The grey shaded regions are identified as - or suspected of - containing metal lines, and are 
thus excluded from the analysis (see Section 4.1 for further details). 
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Table 1. The seven quasar spectra analysed in this study. The columns list, from left to right, the quasar name, redshift and absolute AB 
magnitude at 1450 A, the extent of the region (in kms~^ and proper Mpc) around the quasar where Lya Voigt profiles are fitted to the 
spectrum, and the mean transmission and average signal-to-noise within this region. Further details of the observations and references 
are provided in the remaining columns. 



Name 


y a 
Zq 


A/1450'' 




Rfit 


{F) 


S/N 


Inst. 


Dates 


^cxp 


Refs. 








[km s^-"-] 


[Mpc] 










[hrs] 




SDSS J1148-f5251 


6.4189 


-27.82 


250-3250 


0.3-4.4 


0.471 


23 


HIRES 


Jan 2005-Feb 2005 


14.2 


1 


SDSS J1030-I-0524 


6.308 


-27.16 


250-3450 


0.3-4.8 


0.499 


24 


HIRES 


Feb 2005 


10.0 


1 


SDSS J1623-I-3112 


6.247 


-26.67 


250-3150 


0.3-4.4 


0.573 


21 


HIRES 


Jun 2005 


12.5 


1 


SDSS J0818-I-1722 


6.02 


-27.40 


250-4105 


0.4-6.0 


0.569 


13 


HIRES 


Feb 2006 


8.3 


2 


SDSS J1306-f0356 


6.016 


-27.19 


250-3400 


0.4-5.0 


0.490 


21 


MIKE 


Feb 2007 


6.7 


3 


SDSS J0002-f 2550 


5.82 


-27.67 


250-3150 


0.4-4.8 


0.586 


41 


HIRES 


Jan 2005-Jul 2008 


14.2 


1,3 


SDSS J0836-I-0054 


5.810 


-27.88 


250-3100 


0.4-4.7 


0.556 


27 


HIRES 


Jan 2005 


12.5 


1 



" Redshift and absolute magnitudes are from Carilli et al. (2010). 
References: 

1 - Becker et al. (2006); 2 - Becker et al. (2007); 3 - Calverley et al. (2011) 
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Figure 3. The Doppler parameter CPDFs obtained from the 
near-zones of the seven quasars analysed in this study. The 
Doppler parameters are obtained by fitting Voigt profiles to the 
Lycf absorption in the quasar near-zones, and exclude fits in re- 
gions suspected of containing metal absorption lines. 



each line indicated by the downward pointing arrows. The 
resulting Doppler parameter CPDFs for each of the seven 
quasars are displayed in Figure 3. 



2.2 Hydrodynamical simulations 

The synthetic Lya absorption spectra used to calibrate the 
Doppler parameter CPDF measurements are constructed us- 
ing high-resolution hydrodynamical simulations combined 
with a line-of-sight Lyman continuum radiative transfer al- 
gorithm. The hydrodynamical simulations were performed 
using the parallel Tree-SPH code GADGET-3, which is an 
updated version of the publicly available code GADGET-2 
(Springel 2005). The simulations have a box size of 10/i~^ 
comoving Mpc and a gas particle mass of 9.2 X 10''/i"^Mq. 
Outputs are obtained from the simulations ai z = 6.42, 



z — 6.28, z = 6.01 and z = 5.82. Our procedure for con- 
structing simulated quasar spectra along with the appropri- 
ate numerical convergence tests are discussed in detail by 
BIO. For brevity, we focus instead on describing the addi- 
tional improvements made for this study. 

In order to explore the effect of a wide range of gas 
temperatures on the Doppler parameter CPDF, we perform 
an expanded set of hydrodynamical simulations which em- 
ploy a variety of different thermal histories. We use 14 hy- 
drodynamical simulations in total, which are summarised in 
Table 2. The initial gas temperature in these simulations, 
prior to any quasar heating, is set using a pre-computed 
ultra-violet (UV) background model. The fiducial UV back- 
ground in this study is the galaxies and quasars emission 
model of Haardt & Madau (2001), in which the IGM is 
reionised instantaneously at z = 9. In models A through to 
M, we construct simulations with different initial tempera- 
tures by rescaling the Haardt & Madau (2001) photo-heating 
rates with a constant factor, ^, such that ei — ("ep^"'^. Here 
e™°^ are the Haardt & Madau (2001) photo-heating rates 
for species i = [Hi, Hei, Hen]. 

The temperature at mean density as a function of red- 
shift in a selection of these hydrodynamical simulations, 
prior to any quasar heating, is displayed in Figure 4. The 
photo-heating is coupled to the hydrodynamical response of 
the gas in the simulations; different heating histories will 
therefore result in a different pressure smoothing scale for 
each model (e.g. Gnedin & Hui 1998; Pawlik et al. 2009). 
Since the Lya line widths are sensitive to the changes in 
the gas temperature through Jeans smoothing in addition 
to thermal broadening (Peeples et al. 2010; Becker et al. 
2011), we also consider two further customised UV back- 
ground models, X and Y, for which photo-heating begins at 
z = 12 and z = 15, respectively. These models are displayed 
as the dotted and dashed curves in Figure 4. These extended 
heating histories will result in gas being pressure smoothed 
on larger scales, even if the instantaneous gas temperature 
is similar to our fiducial models. The models are chosen to 
provide a plausible upper limit to the Jeans smoothing scale 
at 2 ~ 6, and we shall use them to explore the effect of 
uncertainties in the IGM thermal history on our results in 
Section 4.3. 
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Table 2. The hydrodynamical simulations used in this work. 
From left to right, the columns list the simulation name, the scal- 
ing factor for the UVB photo-heating rates (see main text for 
details), the redshift which the gas is first photo- heated in the 
hydrodyamical simulations and the volume weighted gas temper- 
ature at mean density, both prior to and after He II photo-heating 
by the quasar. The temperatures in this instance are shown for 
the simulations of SDSS J0818+1722 at z = 6.02, but are broadly 
similar for the other six quasars. 



Model 


c 


^hcat 


log{T^"'"='VK) 


log(To/K) 


A 


0.1 


9.0 


3.29 


4.03 


B 


0.3 


9.0 


3.62 


4.11 


C 


0.5 


9.0 


3.78 


4.16 


D 


0.8 


9.0 


3.93 


4.22 


E 


1.1 


9.0 


4.02 


4.27 


F 


1.8 


9.0 


4.17 


4.35 


G 


2.6 


9.0 


4.28 


4.42 


H 


3.6 


9.0 


4.37 


4.49 


J 


4.7 


9.0 


4.45 


4.54 


K 


5.9 


9.0 


4.51 


4.59 


L 


7.8 


9.0 


4.59 


4.65 


M 


9.3 


9.0 


4.63 


4.69 


X 


Varied 


12.0 


4.00 


4.26 


Y 


Varied 


15.0 


3.92 


4.22 



2.3 Lyman continuum radiative transfer 

We next include the impact of photo-heating by the quasar 
on the surrounding IGM using our line-of-sight radiative 
transfer algorithm. For each of the seven observed quasars, 
we first generate a set of synthetic lines-of-sight using the 
simulation output closest to the quasar systemic redshift. 
A total of 100 lines-of-sight of length 55/i~^ comoving Mpc 
are extracted around the most massive haloes identified from 
each of the hydrodynamical simulations and for each quasar, 
yielding skewers through the IGM density, temperature and 
peculiar velocity field. For our fiducial thermal history, this 
yields a total of 1200 simulated lines-of-sight for calibrating 
each of the seven observed Doppler parameter CPDFs. 

The next stage is to compute the transfer of ionising 
radiation along the lines-of-sight. We assume the quasar 
spectra are described by a broken power law, oc i/""-^ 
for 1050 A < A < 1450 A and oc /z""" for A < 1050 A. 
We adopt a fiducial extreme-UV (EUV) spectral index of 
Qq = 1.5 in this work, consistent with radio quiet quasars 
at lower redshift (e.g. Telfer et al. 2002). The optically 
bright phase of the quasars is assumed to be tq — 10^ yr 
(Haehnelt et al. 1998; Croton 2009). We assume the Hi and 
He: around the z ~ 6 quasars is already highly ionised 
when the quasar turns on [e.g. Wyithe et al. 2008, BIO), 
but that helium has yet to be doubly ionised. The subse- 
quent reionisation and photo- heating of He ii by the quasar 
thus results in an additional IGM temperature increase of 
7000-9000 K. Note, however, an EUV spectral index which 
is softer (harder) than tjq = 1.5 will decrease (increase) the 
amount of Hen photo-heating (Bolton et al. 2009; McQuinn 
et al. 2009). The IGM temperatures at mean density prior 
to and after Hen heating by the quasar are summarised in 
Table 2. 
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Figure 4. The temperature at mean density and its dependence 
on redshift for a selection of the hydrodynamical simulations 
listed in Table 2. The solid curves display the simulations us- 
ing our fiducial thermal history which is based on the galaxies 
and quasars UV background of HMOl. We also explore the ef- 
fect of a more extended period of heating on our results with 
two additional models, X and Y. These are shown by the dotted 
and dashed curves, where the IGM is first heated at 2 = 12 and 
2 = 15, respectively. 



3 METHODOLOGY 

We now briefly review our method for obtaining the near- 
zone temperature constraints from the Doppler parameter 
CPDF (but see BIO for further details). The advantage of 
using the Doppler parameter CPDF is that it fully uses the 
limited number of absorption lines available in the obser- 
vational data and avoids the loss of information associated 
with binning. Although not all the absorption lines in the 
CPDF are purely thermally broadened (Theuns et al. 2000), 
the thermal broadening kernel nevertheless acts on all the 
Lya absorption and the entire CPDF therefore remains sen- 
sitive to changes in the gas temperature. The second ad- 
vantage is that we may compare the observed CPDFs to 
our simulations using the "D-statistic" , which is very sim- 
ilar to the parameter used in a Kolmogorov-Smirnov test 
(e.g. Press et al. 1992). This approach has the advantage 
of providing a non-parametric measure which quantifies the 
difference between Doppler parameter CPDFs drawn from 
different models and the data. On the other hand, obtain- 
ing absolute confidence intervals on the temperature using 
this method requires a large set of realistic simulations for 
calibrating the D-statistic, and these simulations are time 
consuming to perform and analyse. 

For each quasar, we therefore use our sets of synthetic 
spectra to construct Monte Carlo realisations of the Doppler 
parameter CPDF for each of the fiducial simulations, A 
through to M. The simulated CPDFs are obtained by per- 
forming exactly the same line fitting procedure used on the 
observational data. For each simulated line-of-sight in the 
set of 100 spectra, we then measure the D-statistic, which 
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is the maximum difference between the Dopplcr parameter 
CPDF for that simulated line-of-sight and the Doppler pa- 
rameter CPDF for all 100 simulated lines-of-sight: 

Di = max|P(< 6)i - P(< 6)aii|, i = 1...100, (1) 

where we preserve the sign of the difference. The D-statistic 
CPDF for a model with known temperature log To, P(< 
D|logTo), can then be constructed. We then smooth over 
the noise arising from the discrete sampling of the D-statistic 
CPDF with a Gaussian filter of width a = 0.025 before com- 
puting its derivative, ''-P(<-Potgl '"s^o) ^ ]X[q^q that in this work 
we fit 100 spectra for each model, as opposed to only 30 in 
BIO, enabling a finer sampling of the D-statistic distribu- 
tion. 

By then measuring the D-statistic for the observed line- 
of-sight, 7?obs, wc may then use the D-statistic CPDF along 
with Bayes' theorem to obtain a probability distribution for 
the logarithm of the observed temperature at mean density, 
log To: 

p(logToPo..) = ir^^(l^gMZ£)p(iogTo), (2) 

where p(logTo) is the prior on log To and K is a constant 
which normalises the total probability to unity. As in BIO, 
we adopt a flat prior, p(logTo), but adopt a more extended 
prior temperature range, 3.2 < log(ro/K) < 4.8 using our 
expanded set of simulations. The prior is chosen to encom- 
pass the range of full range of initial and final temperatures 
in our simulations (sec Table 2) and is intended to repre- 
sent a reasonable range for the IGM temperature follow- 
ing Hen and/or Hi reionisation (sec e.g. Trac ct al. 2008; 
McQuinn et al. 2009). We then evaluate '"'^^p' at 

twelve discrete points using each of our fiducial hydrody- 
namical simulations, and use a cubic spline to interpolate 
between these to obtain a continuous distribution. Lastly, 
we infer confidence intervals around the median log To by 
integrating p(logTo|T'obB) over the appropriate limits. All 
our temperature measurements are therefore quoted as the 
median of the p(log To | -Dobs) distribution with the confidence 
intervals around the median. This will be implicit through- 
out the remainder of the paper. 



4 SYSTEMATIC UNCERTAINTIES 

We now turn to examine the key systematic uncertainties 
in our analysis. In BIO, we considered five potential sources 
of systematic error. The first three - the mean transmit- 
ted fraction in the quasar near-zone, the influence of strong 
galactic winds, and the dependence of the IGM temperature 
on gas density, were found to have a negligible effect on the 
Doppler parameter CPDF. In the latter case, this is because 
transmission in the 2; ~ 6 near-zones is sensitive to gas close 
to mean density. On the other hand, BIO estimated that 
metal line contamination and uncertainties in the contin- 
uum placement could systematically bias temperature mea- 
surements by up to ~ 2000 K. Although these uncertainties 
arc small compared to the statistical uncertainty on the BIO 
measurement for a single line-of-sight, they will be impor- 
tant to consider for the larger data set used in this study. 
We furthermore examine an additional and important sys- 



tematic which was not included in BIO - the effect of the 
uncertain thermal history at z > 6 {e.g. Becker et al. 2011). 

4.1 Metal lines 

We first examine the effect of metal line contamination on 
our results. As our synthetic spectra do not include absorp- 
tion from intervening metals, wc rrnist remove any metal 
contamination present in the observational data in order to 
avoid biasing our temperature measurements. The erroneous 
identification of metal lines as H i Lya absorption can lead 
to underestimated temperatures, since in general metal ions 
exhibit significantly narrower line widths compared to the 
Lya lines. 

Our approach is to identify possible metal lines in 
the quasar near-zones and exclude these regions from our 
Doppler parameter analysis. Potential contaminants lying 
within the near-zones were first identified by searching for 
metal line systems at other wavelengths in the quasar spec- 
tra, and flagging any associated lines within the near-zones. 
This was followed by flagging any absorption features in the 
near-zones that subsequently remained unidentified and ap- 
peared too narrow to be Hi Lya lines. This procedure is 
most reliable for quasars which have near-IR spectra, en- 
abling greater coverage redward of the Lya emission line. 
Near-IR data was available for five of the seven quasars 
analysed in this work; the exceptions are J1623-I-3112 and 
J1306+0356. 

The contaminated regions identified in the near-zones 
are marked by the grey shaded regions in Figure 2, and 
are excluded from our Doppler parameter CPDF analy- 
sis. For three of the quasars, J0818+1722, J1306-I-0356 and 
J0002+2550, no metal lines were identified. Removing the 
metal contamination slightly raises the temperature con- 
straints obtained from the other four spectra, producing an 
increase of ~ 0.02 dex in our log To constraints. Note, how- 
ever, that it is possible that metal lines which are highly 
blended with the Lya lines remain unidentified. Through- 
out the rest of the paper, we therefore conservatively add 
an additional scatter of 0.02 dex in quadrature to our final 
log To constraints for each line-of-sight. 

4.2 Continuum placement 

The second important systematic identified in BIO was the 
effect of continuum placement uncertainties on the tempera- 
ture measurements. If the continuum is placed too low (high) 
on the observational data, the IGM temperature around the 
quasars will be underestimated (overestimated) as the ab- 
sorption line widths are effectively narrowed (broadened). 
BIO attempted to account for continuum placement uncer- 
tainties by normalising the synthetic spectra by the highest 
transmitted flux in segments of length 1000 km s~^ . The mo- 
tivation for this was to minimise bias by treating the syn- 
thetic data in a similar fashion to the observational data. 

However, the drawback of this approach is that it only 
crudely mimics the continuum fitting process on the ob- 
servational data. In this study, we instead perform "blind" 
normalisations on the synthetic spectra to estimate the con- 
tinuum placement uncertainty (see also e.g. Tytler et al. 
2004; Faucher-Giguere et al. 2008). Our procedure is as fol- 
lows. One of us (JSB) constructed 20 different simulated 
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Figure 5. Left: The continuum level, Cflt , recovered from a blind analysis of 20 synthetic spectra relative to the true unabsorbed continuum, 
Ctrue, against Hubble velocity blueward of rest- frame Lya. The light grey curves display the relative continuum uncertainty for each 
spectrum, while the thick black line shows the average. Upper right: The distribution of the difference between the temperature recovered 
using the algorithm detailed in Section 3 and the input value, logTp"* — logTg" (black solid lines), for 100 synthetic spectra drawn from 
model E where the true continuum is known exactly. The resulting distribution may be approximated by a Gaussian with mean e = —0.01 
dex and standard deviation a = 0.04 dex (red solid curve), indicating that the true temperature is recovered well with some intrinsic 
scatter. The vertical black dotted line and blue dashed lines mark the zero offset and mean of the best fit Gaussian, respectively. Lower 
right: The distribution of log Tq — log Tq" for the same spectra after randomly dividing each by one of the error functions in the left-hand 
panel. The distribution may be approximated as a Gaussian with an offset of e = —0.02 dex and standard deviation a = 0.05 dex. 



spectra with sufficient wavelength coverage on either side 
of the Lya emission line for normalisation. Each spectrum 
was then multiplied with one of the continuum fits compiled 
by Kramer & Haiman (2009), obtained from a large sam- 
ple of low redshift, unobscured quasars. The spectra were 
then processed to resemble the resolution and noise proper- 
ties of the quasars analysed in this work. A second author 
(GDB), who was responsible for continuum fitting the ob- 
servational data, then proceeded to normalise the synthetic 
spectra without any prior knowledge of the true continuum. 

The results of the blind continuum fits to the synthetic 
data are displayed in the left panel of Figure 5, where the 
continuum level. Cat, recovered from the analysis is shown 
relative to the true continuum, Ctruo, against Hubble veloc- 
ity blueward of rest frame Lya. The hght grey curves display 
the relative continuum uncertainty for each of the 20 syn- 
thetic spectra, while the thick black line shows the average. 
On average the continuum is placed to within 5 per cent of 
the true value within ~ 3000 kms~^ of the quasar redshift, 
where the uncertain shape of the quasar Ly-a emission line is 
the largest uncertainty. However, the continuum placement 
is almost always biased low beyond > 3000 km s^^ - by as 
much as 15 per cent on average - where the transmitted fiux 
rarely recovers to the unabsorbed level. Fortunately, since 
the majority of the absorption lines fitted to the observa- 
tional data in this work lie within 3000 kms~^ of the quasar 
systemic redshift, this suggests that the impact of the con- 
tinuum placement on our temperature measurements should 
be relatively modest. 

We quantify this in the right-hand panel of Figure 5, 
where the recovery of the IGM temperature at mean den- 



sity from the simulated Doppler parameter CPDF is tested. 
The upper panel displays the distribution of the difference 
between the input temperature, logTg", and the recovered 
temperature, logTo"*, for 100 simulated spectra drawn from 
model E. This distribution may be approximated by a Gaus- 
sian, displayed by the red curve. In this instance the con- 
tinuum is assumed to be known exactly, and the input tem- 
perature is recovered accurately (to within 0.01 dex, with a 
small scatter of ^ 0.04 dex). 

For comparison, in the lower right hand panel of Fig- 
ure 5 we display the distribution of log Tq — log Tq" for the 
same spectra, but after they have been divided at random 
by one of the error functions displayed in the left hand panel 
of Figure 5. It is apparent that the temperature recovery is 
only very mildly impacted by the continuum placement. The 
recovered temperatures are biased to slightly lower temper- 
atures (by 0.02 dex) with some additional scatter. This is 
a smaller effect than that predicted by BIO, who estimated 
the continuum fitting process could bias the temperature 
measurements by as much as ~ 2000 K. In order to test this 
further, we therefore also performed the same temperature 
recovery test on a set of spectra drawn from model E, but 
this time using the normalisation procedure used by BIO. 
This resulted in temperatures which were biased low by 0.03 
dex. This suggests that the approximate method used in BIO 
slightly over-compensates for the relatively modest system- 
atic offset in the recovered temperature due to continuum 
placement. 

We therefore conclude that the uncertain continuum 
placement will introduce only a small additional uncertainty 
to our results. Nevertheless, we shift all our temperature 
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Figure 6. The effect of the uncertain IGM thermal history at 
2 > 6 on the temperatures recovered from quasar near-zones. 
Top: The distribution of the offset in the recovered temperature 
from the input value, logT^"' — logTp" (black solid lines), for 100 
synthetic spectra drawn from Model X, in which photo-heating of 
the IGM begins at 2 = 12. Bottom: The distribution for model Y, 
with heating beginning at 2 = 15. In both panels, the red curves 
show the best-fit Gaussian, with mean e and standard deviation 
a indicated on the panels, while the vertical black dotted line and 
blue dashed lines mark the zero offset and the mean of the best-fit 
Gaussian, respectively. 



measurements upward by 4-0.02 dex to account for the es- 
timated continuum bias, and a further 0.03 dex is added in 
quadrature to our estimated uncertainties for each line-of- 
sight. 

4.3 Jeans smoothing and the thermal history 

The third and final systematic we examine is the uncertain 
thermal history of the IGM at 2 > 6. This will impact on 
the small-scale structure of the Lya forest at 2: ~ 6 through 
the effect of Jeans smoothing (e.g. Gnedin & Hui 1998; Paw- 
lik et al. 2009). The finite amount of time required for gas 
pressure to respond to changes in the temperature means 
that two models with different thermal histories will be pres- 
sure smoothed on different scales, even if the instantaneous 
temperature is similar at the quasar redshift. In practice, be- 
cause we rely on the simulations to calibrate the relationship 
between the Doppler parameter CPDF and the IGM tem- 
perature, if the true thermal history of the IGM is different 
to that assumed in our models, a systematic bias will be im- 
parted to our measurements. Unfortunately, since the IGM 
temperature is unconstrained at z > 6, the thermal history 
introduces an important uncertainty into our analysis {e.g. 
Becker et al. 2011). 

We investigate the impact of the uncertain thermal his- 
tory on our results in Figure 6, where we examine the ac- 
curacy with which temperatures are recovered from spectra 
drawn from models X and Y. These models have thermal 
histories which are more extended than our fiducial mod- 



els, with heating beginning at z = 12 and z = 15, respec- 
tively (see Figure 4). The increased amount of early heating 
in these two models means that the IGM has had longer 
to dynamically respond to the increased temperatures. As 
a result, the gas distribution in these models is physically 
smoothed on a larger scale. In practice, this means that rel- 
ative to our fiducial simulations, a greater proportion of the 
line broadening in spectra constructed from these simula- 
tions will be due the increased physical extent of the ab- 
sorbers rather than Doppler broadening. 

As a consequence, the recovered distributions for 
log To — log To" in Figure 6 demonstrate that the inferred 
temperatures for models X and Y are significantly higher, 
by 0.07-0.09 dex, compared to the true temperature in 
the models. The measurements also exhibit more scatter 
(cr ~ 0.07-0.09 dex) around the mean than observed for 
the fiducial models (e.g. the upper right panel in Figure 5). 
As might be expected, the systematic offset in the recovered 
temperature is slightly larger for the model where heating 
begins a.t z — 15. The very uncertain thermal history z > 6 
therefore means it is possible that we may overestimate the 
IGM temperature in the real quasar near-zones by as much 
as 0.1 dex. On the other hand, if reionisation occurred very 
late we may underestimate the temperatures. However, the 
latter case will make our conclusions regarding Hen photo- 
heating around the quasars stronger, rather than weaker. We 
therefore address this issue in the next section by presenting 
temperature measurements which include an estimate of the 
impact of additional Jeans smoothing along with our fidu- 
cial results. These Jeans smoothing related uncertainties are 
estimated based on the results of this analysis: we correct for 
a systematic offset by shifting our constraints by —0.08 dex 
and adding an additional scatter of ct = 0.07 dex in quadra- 
ture to the measurement uncertainties for each line-of-sight. 



5 RESULTS 

5.1 Line-of-sight temperature constraints 

We now present the main results of this study by first con- 
sidering the temperature measurements obtained for the in- 
dividual lines-of-sight. The temperature measurements at 
mean density derived for our fiducial thermal history are 
summarised in Table 3, and are plotted as a function of red- 
shift (upper panel) and quasar absolute magnitude (lower 
panel) in Figure 7. It is reassuring to note that the tem- 
perature constraints mirror the Doppler parameter CPDFs 
displayed in Figure 3, with the lowest temperatures de- 
rived from the CPDFs with the smallest Doppler parameters 
on average. The two highest redshift quasars, J 11484-5251 
and J1030-I-0524, exhibit slightly lower temperatures com- 
pared to the rest of the sample, while on the other hand 
J00024-2550 has a slightly higher temperature. However, the 
measurement uncertainties are such that all seven quasars 
are consistent with the same IGM temperature within the 
95 per cent confidence intervals. 

In the upper panel of Figure 7, the temperature 
constraints are also compared to the measurement for 
J0818-I-1722 presented by BIO, shown by the grey diamond. 
Although our revised measurement for J0818-I-1722 is for- 
mally consistent within the 68 per cent confidence interval, 
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Figure 7. Upper panel: The filled circles display the recovered 
near-zone temperatures at mean density against redshift for the 
seven quasar spectra analysed in this work. The thick (thin) error 
bars show the 68 (95) confidence intervals. The data points for 
J1306-I-0356 and J0836+0054 have been offset by Az = -0.03 
for clarity. The grey diamond compares the measurement for 
J0818-I-1722 obtained by BIO. Lower panel: The near-zone tem- 
perature at mean density against quasar absolute magnitude at 
1450 A. 



it is systematically lower by ~ 0.13 dex. The first reason 
for this difference is that in BIO the redshift of J0818-f 1722 
was given as z = 6.00 (Fan et al. 2006). The revised redshift 
oi z = 6.02 presented by Carilli et al. (2010) now extends 
the region where we fit absorption lines by 855kms~^, re- 
sulting in an additional five absorption lines in the near-zone 
(yielding a total of 30). These lines lower the median Doppler 
parameter by ~ 1.7kms~^, and shift the temperature con- 
straint downward by ~ 0.03 dex. The second difference is 
our improved treatment of the continuum uncertainty. As 
discussed in Section 4.2, the approximate continuum uncer- 
tainty correction used in BIO produces temperatures which 
are ~ 0.01 dex higher compared to our new results. The 
third difference is the extended prior probability used in 
this work, especially at lower temperatures. We may ap- 
proximate the BIO measurements by using only models C 
to J for our analysis and restricting the prior probability to 
4.13 < log(ro/K) < 4.56. This yields a temperature ~ 0.03 
dex higher than the value obtained using the full simulation 
set. 



5.2 Average temperature constraints 

The primary improvement of this study is the enlarged sam- 
ple of seven quasars at z ~ 6, which enables us to signifi- 
cantly improve upon the large statistical uncertainty in the 
individual line-of-sight measurement presented by BIO. The 
average temperature constraints derived from the Doppler 
parameter CPDF for all seven quasars are therefore given 
in the final two rows of Table 3. As discussed in Section 4.3, 
each column gives the measurements for both our fiducial 
thermal history and the case of additional Jeans smoothing 
arising from a more extended heating history. 

The average measurements for all seven quasar are dis- 
played in Figure 8, where they are also compared to the IGM 
temperature measurements at z < 5 presented recently by 
Becker et al. (2011). The left hand panel displays the mea- 
surements obtained using our fiducial thermal history, while 
the right hand panel shows our constraints after accounting 
for additional Jeans smoothing. Note that the Becker et al. 
(2011) measurements in the right panel have been shifted 
downward by 2000 K, reflecting the approximate expected 
systematic shift these authors find at z ~ 4-5 when includ- 
ing additional Jeans smoothing in their analysis. 

In each panel, we furthermore show two different ver- 
sions of our 2: ~ 6 temperature measurements. The blue cir- 
cle shows the temperature constraint obtained directly from 
the Doppler parameter CPDF for all seven lines-of-sight. In 
addition, the red square shows the constraint on the ini- 
tial temperature of the IGM before the quasars reionise the 
Hen in their vicinity. Here we have used the Hen heating 
estimates from our radiative transfer simulations to remove 
the contribution of the in-situ Hen photo- heating by the 
quasars. In practice, this is achieved by evaluating Eq. (2) 
using the initial IGM temperature, jjj"'''^''^ in our mod- 
els rather than the temperature after Hen heating by the 
quasar. Note that this constraint assumes the average quasar 
EUV spectral index is Qq = 1.5; a harder (softer) spectral 
index would lower (raise) these constraints. The advantage 
of these measurements is that they attempt to remove the 
impact of Hen heating by the quasar, and should therefore 
more closely represent the temperature of the general IGM 
prior to Hen reionisation. These constraints may therefore 
be more easily compared to expectations for the thermal 
state of the IGM following H i reionisation (Furlanetto & 
Oh 2009; Gen et al. 2009). 

We qualitatively explore the implications of these mea- 
surements for reionisation and the IGM thermal history 
using two simple evolutionary models for the IGM tem- 
perature, shown by the dashed and dotted curves in Fig- 
ure 8. The dashed line, which scales as To oc (1 -I- z)^, is 
normalised to match our temperature constraint including 
Hen heating. This curve represents the maximum possible 
(adiabatic) cooling rate toward lower redshift in the ab- 
sence of any additional heating. In contrast, the dotted curve 
scales as To oc (1 + zf'^^ and is normalised to pass through 
the constraint with Hen heating subtracted. This curve ap- 
proximates the minimum amount of cooling expected if the 
IGM temperature is asymptotically approaching the ther- 
mal state set by the spectral shape of the UV background 
(Hui & Haiman 2003). 

Although the absolute values for the temperatures are 
lower when including additional Jeans smoothing, the gen- 
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Figure 8. Left: The filled blue circle shows the temperature at mean density recovered from all seven quasar spectra analysed in this 
work. The thick (thin) error bars display the the 68 (95) per cent confidence intervals. The red square shows the temperature constraint 
after subtracting the expected photo-heating from the reionisation of Hell by the quasars. For comparison, the open diamonds give the 
temperature measurements from the general Ly« forest at 2 < 5 obtained by Becker et al. (2011) with 2a uncertainties. The slope for the 
IGM temperature-density relation, T = To(l -I- S)"'~^, is assumed to be 7 = 1.3. The dashed curve shows the expected redshift evolution 
if the IGM cools adiabatically, Tq oc (1 -|- z)^, while the dotted curve follows the thermal asymptote, Tq (x (1 -|- z)"'^^. These curves 
approximate the maximum and minimum rate of cooling expected toward lower redshift. Right: As for the left panel, except now the 
near-zone temperature constraints are adjusted for the possible systematic due to additional Jeans smoothing (see Section 4.3 for details). 
Following Becker et al. (2011), the z < 5 measurements have also been shifted downward by 2000 K, refiecting the expected systematic 
shift these authors find at z ~ 4-5 when including additional Jeans smoothing in their analysis. 



eral conclusions we draw from the relative values of the tem- 
perature data remain the same regardless of the uncertain 
thermal history ai z > 6. The results in Figure 8 imply that 
(i) there is an additional and significant source of heating 
around the 2 ~ 6 quasars which is not yet present in the gen- 
eral IGM at z ~ 5 and (ii) there is evidence for a constant or 
gradually increasing temperature in the general IGM from 
z ~ 6 to z ~ 5. The first point suggests that we may be ob- 
serving evidence for the earliest stages of Hen reionisation 
in the immediate environment of these high redshift quasars. 
The second point agrees with the suggestion by Becker et al. 
(2011) that there is a gradual increase in the temperature of 
the general IGM at z < 5 due to the impact of an extended 
epoch of He 11 reionisation globally (but see also the sugges- 
tion by Chang et al. 2011 that significant additional heating 
at z < 5 may arise from TeV blazars, assuming the kinetic 
energy of ultra-relativistic pairs is converted to thermal en- 
ergy via plasma beam instabilities). 

5.3 The thermal proximity effect 

Finally, as we suspect these quasars may have recently 
reionised the Hen in their vicinity, it is interesting to also 
examine the IGM temperature as a function of distance from 
the quasars. We therefore now consider the possibility of di- 
rectly observing the line-of-sight "thermal proximity effect" , 
which will arise from the elevated temperatures expected in 
the Hem regions created by these quasars (Miralda-Escude 
& Rees 1994; Theuns et al. 2002b; Meiksin et al. 2010). 

Assuming the optically bright phase of a quasar is sig- 



nificantly shorter than the recombination timescale, tq <^ 
^rcc^^\ the typical size of a quasar Hem region in proper 
Mpc is 

«„.....Mp.(||)*73i,)'"(l±i)--,3, 

Here N^c — Nucii/W^'' s^^ is the number of Hen ionising 
photons emitted per second, which for the quasars analysed 
in this work ranges from — 1.7-5.3 assuming an EUV 
spectral index of Qq — 1.5. A thermal proximity effect will 
thus be detectable only if the quasar Hem regions are on 
average smaller than the scales over which we make our 
temperature measurements (~ 5 proper Mpc). Note also 
that the adiabatic cooling timescale, which is the appropri- 
ate cooling rate for gas at mean density expanding with the 
Hubble flow, is W ^ l/2H{z) = 7.2 x 10* yr at z = 6 
(cf. the age of the Universe at 2 = 6, fago — 9.6 x 10* yr). 
Consequently, the thermal proximity effect should provide 
a reasonable constraint on the total duration of the opti- 
cally bright lifetime for these quasars. We may therefore 
infer from Eq. (3) that the detection of significantly lower 
temperatures at the edge of the quasar near-zones would 
imply a rather short optically bright phase on average for 
these objects, tq < lO'^ yr. 

We obtain the radial temperature measurements by 
splitting the absorption line fits from all seven quasars into 
four bins of width 750kms~^ over the range 250kms~^ < 
vh < 3250 km s~^. We then use the Doppler parameter 
CPDF in each of the four bins to obtain temperature mea- 
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Table 3. The temperature measurements obtained from the Doppler parameter CPDF in the near-zones of the seven quasars analysed 
in this study. From left to right, each column lists the quasar name, the redshift range over which the Doppler parameters were measured 
from the spectrum, and the final constraints for our fiducial thermal history and for the case of the additional Jeans smoothing expected 
from a very extended period of heating at z > 6 (see Section 4.3 for details). The final two rows also show the line-of-sight averaged 
constraints obtained for all seven quasars, with the last row giving the temperature measurement after subtracting the expected heating 
from the local reionisation of Hen by the quasar. In all instances we assume a flat prior probability of 3.2 < log(ro/K) < 4.8. 
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Figure 9. The temperature at mean density as a function of 
Hubble velocity blueward of the quasar Ly« emission lines. The 
data points are derived from all seven quasars in our sample, 
where the thick (thin) error bars represent the 68 (95) per cent 
confidence intervals. The solid black, dotted red and dashed blue 
curves show the temperature at mean density predicted by ra- 
diative transfer simulations constructed from model D for an op- 
tically bright phase of tq = 10^, 10®-^ and 10® yr respectively. 
Assuming the quasars reionised Hen in their environment, a flat 
thermal profile is inconsistent with an optically bright phase of 
duration tq < 10®-^ yr. 



surements at mean density in the usual manner. Note that 
due to the small number of lines in each bin for individual 
quasars (typically 6 — 8), it is not possible to obtain useful 
constraints from each quasar individually. 

The results of this procedure are displayed as the 
data points with error bars in Figure 9, and are consis- 
tent with constant temperature of log(ro/K) ^ 4.2 within 
^ 3250 kms^^ of these quasars. These are compared to 
a selection of radial temperature profiles from our line-of- 
sight radiative transfer calculations. The simulations are 
constructed using model D, with an initial temperature of 
log(ro/K) ^ 3.9, and are identical to the models used to 
construct our synthetic spectra with the exception that we 
also now consider two shorter optically bright phases of 
tq — lO'''^ yr and tq = lO'' yr. The temperature profiles dis- 
played are averaged over all 100 simulated lines-of-sight and 
smoothed by a box car window of width lOOkms"^ for clar- 
ity. It is evident from this simple comparison that, under 
the assumption the quasars have reionised the Hen in their 
vicinity, the data are inconsistent with an optically bright 
phase with tq < 10®'^ yr. However, due to the relatively re- 
stricted range probed by the near-zones we are unable to 
detect clear evidence for a thermal proximity effect. 

Future progress in this area may be possible at slightly 
lower redshift, where line widths can be analysed at larger 
distances from the quasar. However, detecting the thermal 
proximity effect in quasar spectra at 2 < 3 will likely diffi- 
cult, since it is expected He 11 reionisation is largely complete 
by this time {e.g. ShuU et al. 2010; Worseck et al. 2011). 
Analyses at slightly higher redshift, z > 3.5, may therefore 
be best placed for such a study. On the other hand, a re- 
cent theoretical study by Meiksin et al. (2010) found that 
peculiar motions and IGM density variations can result in 
the Lya line widths exhibiting very little dependence on the 
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distance from the quasar, even in the presence of a thermal 
proximity effect. Detailed simulations which correctly model 
the IGM density and velocity field will therefore be required 
to extract temperature measurements from the data. 



6 CONCLUSIONS 

In this work we present improved measurements of the tem- 
perature of the IGM at mean density around a ~ 6 quasars. 
We use a sample of seven high-resolution quasar spectra, 
combined with detailed simulations of the thermal state of 
the inhomogcneous IGM, to improve upon the first direct 
measurement of the IGM temperature around J0818-I-1722 
at « ~ 6 presented by BIO. This study therefore builds upon 
the work of BIO in three important ways, by using a larger 
sample of quasars, an expanded set of numerical simulations 
for calibrating the temperature measurements, and a more 
detailed analysis of systematic uncertainties. We find that 
the most important systematic is the thermal history of the 
IGM at a > 6, which impacts on our measurements through 
the uncertain contribution of Jeans smoothing to the widths 
of Lyo absorption lines. 

The temperature at mean density averaged over all 
seven lines-of-sight is log(To/K) = 4.21±8:gi (±0.07) at 68 
(95) per cent confidence. On comparison to constraints on 
the temperature of the general IGM at 2: ~ 4.8 which are 
consistent with log(To/K) = 3.9 ±0.1 within 2a, these data 
suggest that there is an additional and significant source 
of heating around the 2 ~ 6 quasars which is not yet 
present in the general IGM at 2: ~ 5. This is most likely 
due to the recent reionisation of He 11 in the vicinity of these 
quasars, which is driven by their hard non-thermal ionising 
spectra. The elevated temperatures may therefore represent 
the first stages of He 11 reionisation in the most biased loca- 
tions in the high rcdshift Universe. It is furthermore found 
that when subtracting the expected amount of He 11 photo- 
heating from the quasars, and assuming a canonical EUV 
spectral index of aq = 1.5, the general IGM temperature is 
similar to that measured at z ~ 4.8. This agrees with the 
suggestion by Becker et al. (2011) that the observed rise in 
the IGM temperature at « < 4.4 is consistent with the onset 
of an extended epoch of He 11 reionisation globally. 

We also examine the evidence for a line-of-sight ther- 
mal proximity effect around these quasars by analysing the 
Doppler parameters for all seven lines-of-sight in radially 
spaced bins. We find no clear evidence for a thermal prox- 
imity effect due to an He iii region around the quasar, but 
note that the limited extent of the near-zone prevents us 
from detecting photo- heated Hem bubbles larger than ~ 5 
proper Mpc in size. Under the assumption that the quasar 
has reionised the He 11 in its vicinity, the data are therefore 
inconsistent with a short optically bright phase tq < l(f '^yr. 

Finally, in this work we have not examined the implica- 
tions of our temperature measurements for Hi reionisation 
at 2: > 6. These high-rcdshift temperature measurements 
should still probe the thermal signature of this landmark 
event, potentially yielding valuable insights into the timing 
and duration of this process {e.g. Theuns et al. 2002a; Hui 
& Haiman 2003; Furlanetto & Oh 2009; Trac et al. 2008; 
Cen et al. 2009). We intend examine this in detail in future 
work. 



ACKNOWLEDGMENTS 

The hydrodynamical simulations used in this work were 
performed using the Darwin Supercomputer of the Uni- 
versity of Cambridge High Performance Computing Service 
(http://www.hpc.cam.ac.uk/), provided by Dell Inc. using 
Strategic Research Infrastructure Funding from the Higher 
Education Funding Council for England. JSB acknowledges 
the support of an ARC Australian postdoctoral fellowship 
(DP0984947), and GDB thanks the KavU foundation for fi- 
nancial support. 



APPENDIX A: TESTS OF THE 
TEMPERATURE MEASUREMENT 
PROCEDURE 

Tests of the temperature measurement procedure described 
in Section 3 are displayed in Figure Al. The solid black lines 
show distributions of the difference between the temperature 
at mean density measured from 100 synthetic spectra and 
the input value, logTo°"* - logT,5°, for four of our fiducial 
models. The models span a temperature range consistent 
with the observational measurements presented in this work. 
The distributions may be approximated by a Gaussian with 
mean e and standard deviation a, shown by the solid red 
curves. The means are within < 0.02 dex of the true input 
temperature in all cases, indicating our temperature mea- 
surement procedure is reliable to this level in the absence 
of any additional systematic uncertainties. Note, however, 
the accuracy of the temperature recovery is slightly poorer 
for model K, where the distribution is skewed to lower tem- 
peratures. This will be true in general for models with pro- 
gressively higher temperatures; the line widths scale as the 
square root of the temperature, b oc T^''^, and the Doppler 
parameter CPDF thus has less discriminative power as the 
temperature increases. 
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